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We propose a computational scheme for materials where standard Local Density Approximation 
(LDA) fails to produce a satisfactory description of excitation energies. The method uses Slater's 
"transition state" approximation and Wannier functions basis set. We define a correction to LDA 
functional in such a way that its variation produces one-electron energies for Wannier functions 
equal to the energies obtained in "transition state" constrained LDA calculations. In the result 
eigenvalues of the proposed functional could be interpreted as excitation energies of the system 
under consideration. The method was applied to MgO, Si, NiO and BaBiOa and gave an improved 
agreement with experimental data of energy gap values comparing with LDA. 

PACS numbers: 78.70.Dm, 71.25. Tn 

I. INTRODUCTION 

The one-electron eigenvalues in LDA calculations are Kohn-Sham energies that were defined as auxiliary Lagrange 
multipliers in the problem of Density Functional minimizationi. As such they formally can not be used to calculate 
spectral properties of the system because they do not directly correspond to excitation energies. However for the cases 
where electron states are itinerant, for example wide band metals, excitation spectra calculated with LDA eigenvalues 
were found to be in good agreement with experimental photoemission and optical spectra. The agreement is much 
worse for band insulators and semiconductors where LDA gives systematically underestimated values of the energy 
gap^. For Mott insulators, for example transition-metal oxides, LDA calculated spectra could be qualitatively wrong, 
giving metallic state while experimentally those systems are wide gap insulators^. 

There were many attempts to cure this deficiency of LDA. Among the most widely used methodes one can mention 
GW— , SICjS and LDA+U§i While those approaches have their advantages there is still no universally accepted 
calculation scheme which would be as simple and practical as standard LDA and a search of better methods continues 
in scientific community. 

The basic problem of a Density Functional Theory (DFT)(and consequently of the LDA as one of it's approxima- 
tions) is that DFT was designed to reproduce a ground state properties of a system. In order to obtain excitation 
energy in the spirit of DFT one must calculate total energy of the system in the excited state. Then its difference from 
the ground state total energy corresponds to excitation energy. The excited state can be calculated minimizing LDA 
functional with the constraining condition that occupancy of specific one-electron state differs from its occupancy in 
the ground state. Such "constrained LDA" calculations indeed could give good results for excitation energies^. 

Slater had shown& that a good approximation for excitation energy can be one-electron energy for the specific 
one-electron state calculated in a so called "transition state", where its occupancy is equal 0.5 being half way between 
the values for the final and initial state of the excitation process. Transition state calculations indeed gave excitation 
energies for free atoms and ions in a good agreement with experimental data&. 

However this approach can not be directly applied to extended systems such as solids with translation symmetry 
where one-electron states arc Bloch functions extended over infinite crystals. The change of their occupancies produces 
infinitesimal change of electron density and so constrained LDA calculations for "transition state" give exactly the 
same onc-clcctron energy values as in a ground state. 

A full set of Bloch functions can be replaced by the equivalent set of Wannier functions produced by some unitary 
transformation (see Eqs. and J2J). In contrast to Bloch states Wannier functions are localized in a specific space 
area. If one choose Wannier functions as one-electron states in "transition state" constrained LDA calculations then 
there will be a significant change in the charge density comparing with a ground state and hence a corresponding 
correction to the excitation energy. 

The idea of our "generalized transition state method" developed in the present work is to use "transition state" 
approximation to define a modified functional in such a way that Wannier functions one-electron energies calculated 
with this functional would have a meaning of excitation energies for electrons on the corresponding Wannier function 
states. This is achieved by adding to LDA functional a correction changing Wannier functions energies from LDA 
values to those obtained in the "transition state" calculations. 

The main effect of " transition state" calculations comparing with standard LDA is lower one-electron energy values 
for the occupied states and higher values for the unoccupied states due to the decreased occupancy of the former and 
increased one for the latter. Then application of the correction should enhance energy separation between occupied 
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and empty bands and hence give an increased value for the energy gap thus improving agreement with the experimental 
data. 

To test our method we chose a band insulator MgO, semiconductor Si, Mott insulator NiO and Peierls insulator 
Ba.Bi03. In all those systems LDA gave strongly underestimated energy gap values. Calculations by "generalized 
transition state method" resulted in band structure for all those materials being in a good agreement with experimental 
spectroscopy data. 

The paper is structured as follows. In Sec.|n]we give the details of our calculation scheme. In Sec. lIIII the results for 
the electronic structure of MgO, Si, NiO and BaBi03 obtained by the method developed in this work are presented 
and compared with the standard LDA calculations and experimental data. Finally in Sec. II VI we close this work with 
a conclusion. 



II. METHOD 

The orbital projection calculation scheme for Wannicr functions (WF's) used in the present work was described in 
details in the earlier paper— where the LDA+DMFT (DMFT - Dynamical Mean-Field Theory) method in Wannier 
function basis set was proposed. Below we present the main formulas of this scheme. 

A. Definition and construction of Wannier functions 

The concept of WF's has a very important place in the electron theory in solids since its first introduction in 1937 
by Wannier WF's are the Fourier transformation of Bloch states l^k) 

i^ T ) = ^E e ~' ikT i^>, (!) 

where N is the number of discrete k points in the first Brillouin zone (or, the number of cells in the crystal) and T 
is lattice translation vector. 

Wannier functions are not uniquely defined for a many-band case because for a certain set of bands any orthogonal 
linear combination of Bloch functions IV'ik) can be used in Eq. Q). In general it means that the freedom of choice of 

Wannier functions corresponds to freedom of choice of a unitary transformation matrix for Bloch functions^ 

i&k) = £e#W>- ( 2 ) 
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There is no rigorous way to define f/jj • This calls for an additional restriction on the properties of WF's. Among 
others Marzari and Vanderbiltii used the condition of maximum localization for WF's, resulting in a variational 
procedure to calculate Uf^ . To get a good initial guess the authors ofii proposed choosing a set of localized trial 
orbitals \(f> n ) and projecting them onto the Bloch functions |"0ik). It was found that this starting guess is usually quite 
good. This fact later led to the simplified calculating scheme ir±i£ where the variational procedure was abandoned 
and the result of the aforementioned projection was considered as the final step. 

For the projection procedure used in the present work one needs to identify the set of bands and corresponding set 
of localized trial orbitals \<j> n ). The choice of bands and orbitals is determined by the physics of the system under 
consideration and will be discussed later. 

The set of bands can be defined either by the band indices of the corresponding Bloch functions (Ni, iVjj), or by 
choosing the energy interval (E\,E2) in which the bands are located. Non-orthogonalizcd WF's in reciprocal space 
|W„k) are then projection of the set of site-centered atomiclikc trial orbitals \<p n ) on Bloch functions |V>ik) of the 
chosen bands defined by band indices (Ni to N%) or by energy interval (Ei,E 2 ): 

N 2 

\Wnk) = E l^k)(^k|0n> = IVtkXV'ik I'M- (3) 

»=JVi i(E 1 <e i (k)<E 2 ) 

In the present work we have used LMTO method^ 3 - to solve a band structure problem and the trial orbitals \4> n ) were 
LMTO's. Note that in the multi-band case a WF in reciprocal space |W„k) does not coincide with the cigenfunction 
|^nk) due to the summation over band index i in Eq. ©. One can consider them as Bloch sums of WF's analogous to 
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the basis function Bloch sums <?!> k (r) (Eq. iJSJ). The coefficients (ipik\(j) n } in Eq. Q define (after orthonormalization) 

the unitary transformation matrix in Eq. (§2J) ■ 

In any DFT method the Kohn-Sham orbitals are expanded through the certain basis: 

= E c M>- ( 4 ) 

The basis functions of the LMTO method are Bloch sums of the cite centered orbitals: 

#M = VX7E e * kT <Mr-R 9 -T), (5) 



where \i is the combined index representing qlm (q is the atomic number in the unit cell, Im are orbital and magnetic 
quantum numbers), R g is the position of atom in the unit cell. 
For the orthogonal LMTO basis = (0 M |V'tk) and hence 

\w nk ) = £ i^ k ) C ^ = ^ E c Mi c »*i^) = E^i#>> ( 6 ) 

with 

b k = c k c k * (7) 

For a nonorthogonal basis set orthogonalization of the Hamiltonian must be done before using Eq. Q) . 
In order to orthonormalizc the WF's Eq. © one needs to calculate the overlap matrix O nn /(k) 

N 2 

O nn ,(k) = (W nk \W n , k ) = £ (8) 

i=JVi 

and its inverse square root S nn > (k) 

s nn ,(k) = o m y 2 (k). (9) 

In the derivation of Eq. (JSJ the orthogonality of Bloch states (VViklVVk) = <5rm' was used. 
From Eqs. © and the orthonormalized WF's in k space |W„k) can be obtained as 

\w nk ) = £ s nn ,(k)\w n , k ) = £ i^k)^.* = J2 *£J#>. ( 10 ) 

n' i—Ni f-i 

with 



iV 2 

6 k „ = (^|W„ k ) = ^ (12) 

i=Ni 

The real space site-centered WF's at the origin \W°) are given by the Fourier transform of |W n k) (Eq. Q) with 
T = 0. From Eqs. (JTUJ) and JSJ one finds 



W„(r) 



iEw^) - E(^E^«W- T ) 

V k T> V k 7 



^ W '(n,/i,T)^(r-T) (13) 



= ^2w(n,s)(l> a{s) (r -T s ), 
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where w' and w are the expansion coefficients of WF in terms of the corresponding LMTO orbitals, in particular. 

^' s ) = ^E eJkT ^( sK (14) 

k 

Here s is an index counting the orbitals of the neighboring cluster for the atom where orbital n is centered (T s is the 
corresponding translation vector, a(s) is a combined qlm index). The explicit form of the real space WF Eq. (|13[) can 
be used to produce, e.g., shapes of chemical bonds. For other applications only the matrix elements of the various 
operators in the basis of WF Eq. I|1(J|) arc needed. 

Using Eqs. fTT))l . ifTTjl. and (fT^)) one can find energies of WF's: 



E™ F = (^„ T |^|^k) £ ,(k)(^ k |j|^) 



TV 

k i=Ni 



and their occupancies: 

QT = {Wj\(Y,\M9(.E f -e t (k))(^ k \)\Wj 



= ^EE ZiftAEj-eiM), (16) 

k i=Ni 

where ej(k) is the eigenvalue for a particular band, 9{x) is the step function, Ef is the Fermi energy. 

The transformation from LMTO to WF basis set is defined by the explicit form of WF's Eqs. 1)10(1 . i|12|) . and by the 
expressions for matrix elements of the Hamiltonian and density matrix operators in WF basis (Eqs. i|15|) an d l|16|ll. 
The transformation from WF to LMTO basis can also be defined using Eq. (|1()J) . Such transformation is needed 
in calculations using correction potential in the form of Eq. (|25[) and for constrained LDA calculations determining 
"transition state" energies of WF's. For example if constrain potential is diagonal in WF basis (H nn i = A n S nn '), then 
in LMTO basis its matrix elements can be calculated via: 

H constr =Y,\ W n) A n(W?\ (17) 



\W 



E 



(18) 



H^(k) = (cf> k jH constr \^) (19) 



E 



b k A b k * 



B. "Generalized transition state" method 



Transition state calculation scheme proposed by Slater— allows to calculate excitation energy for the process of 
adding (removing) an electron to (from) the system from (to) the infinity where potential is supposed to be equal 
zero. For that one should calculate LDA eigenvalue (Kohn-Sham equations eigenvalue) of the corresponding one- 
electron state with its occupancy equal 0.5. In other words for the occupied states the occupancy is reduced by one 
half and for the empty ones it is increased by one half. 

This scheme can be derived in the following way. The total energy difference between final and initial states for 
the process of electron addition to the one-electron state j can be calculated as an integral of total energy derivative 
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over occupancy qj (this derivative ism the corresponding LDA eigenvalue ej — ^ E -) 



1 

f { OE 

E( qj = 1) - E(q 3 = 0) = J dq 3 1^ — 

o J 

1 

dqjCjiqj) w £j (0.5) 



(20) 



n 



It means that the one- electron energy (LDA functional eigenvalue) calculated with its occupancy value half way 
between final and initial states is a good approximation for the total energy difference (excitation energy). This 

equation becomes exact if LDA eigenvalue £j(qj) is a linear function of occupancy qj ( — const) , which is 

usually with a good accuracy true. In this case for empty state (qj = 0): 

!j( 0.5) = ej (0) + i»|M (21) 

For occupied states (<jy = 1) a sign plus in Eq. (|21(l will be replaced by minus. The general formula will be: 

6,(0.5) = ej (q j ) + ( 1 --q j )^l (22) 

So effect of "transition state" correction to LDA values is to increase energy for empty states (addition energy) and to 
decrease it for the occupied states (removal energy). That will also results in a larger value of energy for the excitation 
from occupied to empty states of the system. 

Transition state method can be reformulated in a functional formalism. For that equation (|22J) should be obtained 
via variation of an auxiliary functional. This "transition state" functional Ej<s is defined by adding to LDA functional 
Elda[p) the correction termi^ depending on the occupancies qj: 

E TS = E LDA [ P ] - \ £ ~ \f (23) 

j 3 

All eigenvalues obtained from this functional will automatically have a correction (i — qj) to LDA eigenvalues 

like in the right part of Eq. 11221) . Then one calculation for the ground state of the functional Ets gives result equivalent 
to a set of "transition state" calculations for every one-electron state. However one still needs to run a set of constrained 
LDA calculations in order to determine derivatives ^q^ 3 ^ ■ 

The " transition state" method to calculate excitation energies has proved to be sufficiently successful for small size 
systems like atoms and ions£ but it can not be directly applied to solids. Bloch functions are extended over the crystal 
and a change of their occupancy will have a negligible effect on the charge density. In this case "transition state" 
calculations will give the same one-electron energies as the ground state results. 

For a case of fully occupied bands transformation from a set of Bloch states to Wannier functions does not change 
charge density distribution. It is more than that, many-electron function defined as a Slater determinant constructed 
from one-electron Bloch functions is identical to the Slater determinant made out of Wannier functions. (The de- 
terminant value is not changed in the result of adding to one of its rows a linear combination of other rows and 
Wannier functions are by definition linear combinations of Bloch functions (see Eqs. Q and J2J)))- Then ground state 
properties are invariant to such transformation. 

However for excited states there is an important difference between Bloch and WF representations. As WF is 
localized, a change of its occupancy will result in a finite charge density modification in the area of its localization 
and hence " transition state" calculations will result in a significant correction for the WF one-electron energy. While 
in standard LDA one-electron energies arc Kohn-Sham energies that were defined as auxiliary Lagrange multipliers 
in the problem of minimizing Density Functional, WF energies obtained in "transition state" calculations have a 
meaning of excitation energies for an electron on the corresponding WF states. 

The idea of our "generalized transition state" method (GTS) is to replace in the original formulation of "transition 
state" method (Eqs. 12012211 s ) a set of one-electron states in a form of Bloch functions for an infinite crystal by an 
equivalent set of Wannier functions. The corresponding functional will be analogous to Eq. lO but with occupancies 
Q^ F (Eq. H16|) ~) and energies E^ F (Eq. (|1 5|) ~) corresponding to Wannier functions: 

1 dE WF 1 
E G ts = E LDA - - £ qqwf ^ F ~ 2^ ^ 
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The variation of the correction term in (|24|l will produce a correction Hamiltonian H corr in the form of projection 
operator: 

H corr = Y.\ W nWn{W^\. (25) 
nT 

Then Hamiltonian operator for "generalized transition state" method is given by: 

HGTS = HlDA + H corr . (26) 

\Wn) in Eq. (|25J) are Wannier functions (Eq. I|18(l ) and 6V n are defined as a difference between Wannicr functions 
one-electron energies calculated in "transition state" and ground state (see Eq. 122ft : 

«t-!w<j-«r') pt) 

The values of derivatives d( ^iv f (or equivalently SV n themselves) should be determined in constrained LDA calcula- 
tions. 

The calculation scheme is the following (see Fig^l. To define correction operator (Eq. I|25fl) one needs to know 
a set of "transition state" corrections to WF's energy values SV n (Eq. and explicit form of WF's determined 

by expansion coefficients in the basis orbitals (LMTO) (Eq. flUQ - Both SV n and b^ n should be calculated 
self-consistently. On every self-consistency loop Bloch functions \ipiu) calculated with "generalized transition state" 
method Hamiltonian Hqts (Eq. (|26fl 'l are used to define new Wannier functions |W„k) vi a Eq. © to get a new set 
of coefficients b^ n . Then a scries of constrained LDA calculations (using constrain potential in the form of Eq. (|17|) ) 
for every type of WF | W®} (Eq. $13$ ) is performed where the occupancy (Eq. i|16|n of this particular WF is kept to 
be 0.5. The energy of this WF's is then calculated using Eq. I |15L The new SV n is defined as a difference between 
this "transition state" value and the corresponding value from calculation where the WF occupancy is the same as in 
ground state. 



The important question is what set of bands and trial orbitals should be used in projection procedure (Eq. Q3J) to 
calculate WF's. In band insulators and semiconductors valence band corresponds to bonding states and conduction 
band to antibonding ones. So if one will calculate WF's via projection procedure separately for valence and conduction 
bands (summation over band indices in Eq. (|3J) is running over occupied bands only for valence WF's and over 
empty bands for conduction states WF's) then the results will be bonding and antibonding functions extending over 
neighboring atoms. But if WF's will be obtained using a full set of valence and conduction bands (summation over 
band indices in Eq. is running over combined set of occupied and empty bands both for valence and conduction 
states WF's), then the opposite sign contributions on neighboring atoms from bonding and antibonding functions 
cancel each other and the resulting WF resembles an original atomic orbital. 

For spectral properties where excitation occurs from (to) valence (conduction) states WF should represent the 
corresponding bonding (antibonding) functions. Then the most natural choice would be to use in projection procedure 
(Eq. (jSJ)) two separate sets of occupied and empty bands and two sets of atomic- like orbitals that give a strongest 
contribution to the corresponding bands. For example in MgO case the orbitals for occupied bands would be oxygen 
2p and for empty bands Mg 3s and 3p orbitals. The application of the correction operator (Eq. (J5SJ)) will shift relative 
energies of the bands but will not change significantly wave functions themselves. They will remain the same as were 
obtained in LDA calculations. Then charge and spin electron density will be not modified by GTS correction. 

This approximation could be satisfactory for semiconductors and band insulators where the only problem of LDA 
is underestimation of excitation energies. However for Mott insulators LDA calculations could give even ground state 
properties qualitatively wrong. For example for parent high-T c compound La2Cu04 LDA gives nonmagnetic metallic 
solution while experimentally it is an antifcrromagnctic insulator—. For NiO LDA gave strongly underestimated value 
of magnetic moment on Ni ion (1.0 fiB versus experimental value sal. 8-1. 9 ^b^)- The problem with these systems 
could be traced to underestimation in LDA of the energy of virtual excitations from occupied oxygen 2p to empty 
transition-metal 3d orbitals. This leads to overestimation of hybridization between those orbitals and hence decreases 
a tendency to magnetic moment formation. 

In order to correct the LDA underestimation of virtual excitations energies one must use in "transition state" 
calculations WF's resembling pure atomic orbitals such as oxygen 2p and transition- metal 3d orbitals and not their 
bonding (antibonding) combinations which are valence (conduction state) WF's. As it was explained above for that 
one should choose in projection procedure (Eq. (JjJJ) the combined single set of valence and conduction bands instead 
of the two separate sets of occupied and empty bands. The "transition state" calculations for these Wannier functions 



7 



n start tt 
P , ti 



LD A start 



k start 



8V. 



start 
n 



SVn 



Series of constrained 
calculations for 
each type of WF 



H 



GTS 



TjLDA | TT 



corr 



p,H 



LDA 



[in 



Construction of atomic-like 
WF's using combined set of 
valence and conduction bands 



LDA 



5V n 



Series of constrained 
calculations for 
each type of WF 



Construction of valence and 
conduction states WF's using 
filled and empty bands 



H 



GTS 



H LDA + H c 



Excitation properties 
(bands, DOS] 



FIG. 1: The calculation scheme of "generalized transition state" method. 



will give the result which have a meaning of energies for virtual excitations between atomic orbital states and not real 
excitations which happen in spectroscopy experiments from (to) valence (conduction) band states. The corresponding 
correction in the form of Eq. 125f) will increase the energy of virtual excitations from occupied oxygen 2p to empty 
transition-metal 3d orbitals and hence decrease hybridization between those orbitals that should enhance a tendency 
to magnetic moment formation. 

That means that a different definition of WF's is needed in calculations for ground state properties (atomic orbital 
WF's obtained by projection on the combined single set of valence and conduction bands) and in calculations for 
spectral properties (valence (conduction state) WF's projected by using separately valence bands for occupied states 
and conduction bands for unoccupied ones). 

In the present work we performed two steps of calculation for every system (see Fig^l . In the first one (part of the 
FigUscheme bounded in a bold line rectangle) we chose a single combined set of bands to define WF's resembling pure 
atomic orbitals. Calculations arc repeated till achieving full self-consistency in charge and spin densities p(r) (defining 
one-electron LDA potential and hence parameters of LDA Hamiltonian Hlda)i s °t of "transition state" corrections 
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TABLE I: Comparison of calculated and experimental energy gap values (eV) 

to WF's energies SV n and WF's themselves determined by expansion coefficients b^ n (Eq. (|12fl ). In the result ground 
state properties of the system are calculated which are modified from those in LDA solution via "transition state" 
correction to energy parameters of virtual excitations. The next step is to calculate "transition state" correction to 
the energies of real excitation observed in spectroscopic measurements. At this stage we fix obtained on the first step 
Bloch functions \ipik) and parameters of LDA Hamiltonian Hlda- From these Bloch functions IV'ik) new valence and 
conduction state Wannier functions are calculated (new expansion coefficients b^ n (Eq. dJ))) in projection procedure 
(Eq. ||3J)) using two separate sets of occupied and empty bands. These new WF's define a new set of "transition 
state" corrections to WF's energies SV n via constrained LDA transition state calculations and hence a new correction 
Hamiltonian (Eq. (|25fl ) which is used to calculate spectral properties of the system. 

The above described calculation scheme with two different definitions of Wannier functions could seem to be 
overcomplicated, but we have found that no single set of WF's can describe both ground state and excitation properties. 
In general, the choice of procedure defining Wannier functions is dictated by the physics which we want to describe 
in our calculations. 



III. RESULTS AND DISCUSSION 

In order to test proposed above "generalized transition state" (GTS) method we have performed calculations for 
four systems representing various types of electronic structure: simple metal oxide band insulator MgO, covalent bond 
semiconductor Si, transition metal oxide Mott insulator NiO and Peierls insulator BaBiOa. In all those cases we have 
obtained significant change of electronic structure leading to good agreement of calculated and experimental energy 
gap values (see Table Pi. 

A. MgO 

Simple metal oxide MgO is a good example of band insulator. It has a cubic NaCl crystal structure with lattice 
parameter equal to 4.21 A. Experimental value of energy gap (7.83 eV— ) is nearly 3 eV larger then LDA calculated 
value (5.04 eV in our calculations in good agreement with other worksSi). The choice of energy bands and orbitals 
needed for calculation of Wannier functions via Eq.© is straightforward in this case: oxygen 2p orbitals for three 
occupied valence bands and magnesium 3s and 3p orbitals for unoccupied conduction bands. Constrained LDA 
calculations for transition state gave the following corrections to the energies of Wannier functions: 5Vmg3s = 1-44 
eV, 6Vm 9 3 P = 1.81 eV, SVo2p = —2.41 eV . Stronger effect of GTS correction for 02p states comparing with Mg3s, 3p 
WF's could be understood taking into account more extended nature of magnesium orbitals comparing with oxygen 
ones. On Fig[21 band structure calculated by GTS method is presented together with the results of standard LDA 
calculations. According to the SV correction values valence bands formed by 02p states are pushed down and 
conduction Mg3s, 3p bands arc pushed up increasing the energy gap value to good agreement with experimental data 
(see Table IJ). Please note that GTS correction potential (Eq. tyity ) is not a rigid shift of bands and so the difference 
between LDA and GTS gap values is not simply given by the corresponding difference of SV correction values. 

B. Si 

While MgO is an ionic compound, silicon represents a simplest case of a covalent bond semiconductor. Si has 
diamond crystal structure with lattice parameter equal to 5.43 A. For this system we used four types of Wannier 



9 




- 




L G X W L K G 



FIG. 2: MgO band structure. Dashed lines shows LDA results and solid lines correspond to GTS calculations 

functions: SiSsi and SiSpi calculated for the valence bands and Si3s2 and Si3p2 for conduction bands. The corre- 
sponding values of GTS potential (Eq. (|25p) corrections obtained in constrained LDA transition state calculations 
are: SV Si 3 Sl = -0.82 cV, 5V S a Pl = -0.44 eV, SV Sl 3 S2 = 0.28 eV, 8Vs%a P2 = 0.71 eV. The resulting band structure is 
shown on Fig|3] Negative potential correction for valence bands and a positive one for conduction bands produce a 
sizable increasing (0.60 eV) of the energy gap value resulting in a good agreement with experimental data (sec Table 

13. 

C. NiO 

In contrast to energy gap values the ground state properties for MgO and Si are satisfactory reproduced by LDA 
and the first step of GTS calculations with atomic-orbital-like WF's (see Sec. Ill Bl and Fig0 did not lead to significant 
changes comparing with pure LDA. However for transition metal oxides even the type of the ground state can be 
given wrong by LDA, like for cupratcsi For nickel oxide NiO (cubic NaCl crystal structure with lattice parameter 
equal to 4.17 A) the LDA error is not so severe: the LDA solution is an antifcrromagnetic insulator in agreement with 
experiment. However LDA magnetic moment value for Ni ion is strongly underestimated: 1.0 fiB versus experimental 
value wl. 8-1.9 ^b^- In this case the first step of GTS calculations responsible for the "transition state" correction to 
energy parameters of virtual excitations was essential: calculated magnetic moment value was increased in the results 
of GTS correction from 1.0 [is to 1.8 fiB in good agreement with experimental data. There were five different types of 
Wannier functions in calculations for NiO: Ni3d-^2 3 1\Ni3d-e g j, Ni3d-<2 9 I, Ni3d-e g 1 (the only unoccupied states for 
NiO) and 02p (by symmetry 02p states are not spin-polarized). Constrained LDA calculations for transition state 
gave the following corrections to the energies of Wannier functions:(5VAri3d-t, s t = —1-97 cV, 5VNi3d-e g t = —1-97 eV, 
SV Nm -t 2g l = - 2 -39 eV, SV Ni 3d-e g l = 2.11 cV and SV 02p = -1-14 cV. 

The energy band dispersions obtained in GTS calculations together with LDA bands are shown on FigQ] The energy 
gap value was increased dramatically: from tiny 0.11 eV to a 3.76 eV value in good agreement with experimental data 
(see Table HI). Not only the energies of unoccupied Ni3d — e g j states were pushed up and those for occupied bands 
pushed down as it was the case for MgO and Si. The different values of SV correction for occupied Ni3d and 02p 
states result in a smaller energy separation between the corresponding bands in the occupied part of the calculated 
DOS (sec FiglSJ). Note also decreased strength of admixture of Ni3d — e g J. states to the oxygen band energy area 
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FIG. 3: Si band structure. Dashed lines shows LDA results and solid lines correspond to GTS calculations 

resulting in a more ionic nature of GTS solution comparing with pure LDA and hence magnetic moment value 1.8 
Hb closer to pure ionic value 2.0 /is 

D. BaBi0 3 

BaBi03 is an interesting example of the Peierls insulator or "negative U" system^. Formal valency of bismuth is 
+4 and this corresponds to the half-filled Bi6s shell. Stable valent states for Bi arc +3 and +5 and those chemical 
arguments are often used to explain the experimental distorted cubic perovskitc crystal structure of BaBiC>3^. In 
addition to the tilting of BiOg octahedra there is also so called "breathing" distortion producing inequivalent Bil and 
Bi2 crystallographic positions with expanded and contracted Bi-0 bond lengths. Bil can be associated with Bi +3 
and Bi2 with Bi +5 . This distortion leads to opening an energy gap so BaBiC>3 can be seen as a three-dimensional 
Peierls insulator. The Fermi level in LDA calculations for ideal cubic perovskite crystal structure crosses half-filled 
band of Bi6s symmetry which can be described by an effective half-filled Hubbard model. An experimentally observed 
instability toward formation sites with empty (Bi +5 ) and completely filled (Bi +3 ) Bi6s shell can be interpreted as a 
negative value for a Coulomb interaction parameter U. 

LDA calculations for BaBi03 gave a negative or very small positive values for energy gap^S while spectroscopy 
measurements demonstrate a sizable value of «0.5 eV—. For our GTS calculations we have chosen a minimal set of two 
Wannier functions calculated using orbitals of Bil-6s and Bi2-6s symmetry and bands located in the energy window of 
±5 eV around Fermi energy. Constrained LDA calculations for transition state gave the following corrections to the 
energies of Wannier functions: 5Vbh-6s = —0.19 cV (5Vbj2-6s = 0.43 eV. The energy band dispersions obtained in 
GTS calculations together with LDA bands are shown on Fig|SJ Only two bands close to Fermi energy were effected 
by GTS correction and energy gap value has increased from 0.15 eV to a 0.48 cV value in good agreement with 
experimental data (see Table UJ. 
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FIG. 4: NiO band structure. Dashed lines shows LDA results and solid lines correspond to GTS calculations 



IV. CONCLUSION 

We developed a calculation scheme based on "transition state" idea of Slater and Wannier functions set of one- 
electron states. This method was applied to the four materials representing ionic compounds (MgO), covalent bond 
semiconductors (Si), transition metal oxide Mott insulators (NiO) and Peierls insulator BaBiC>3. The results have 
shown a systematic improvement of energy gap values. Not only excitation energies but also ground state properties 
such as magnetic moment value for NiO can be significantly improved comparing with LDA results. Encouraged by 
the promising results reported in the present paper we plan to apply this method to other materials where standard 
LDA approach fails. Such calculations are in progress. 
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